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CN Abstract 

^ We study the fully discrete elliptic integrable model Q4 and its immediate trigonomet- 

ric and rational counterparts (Q3, Q2 and Ql). Singular boundary problems for these 
equations are systematised in the framework of global singularity analysis. We introduce a 
technique to obtain solutions of such problems, in particular constructing the exact solution 
on a regular singularity-bounded strip. The solution technique is based on the multidimen- 

»— ( sional consistency and uses new insights into these equations related to the singularity 

structure in multidimensions and the identification of an associated tau-function. The 
obtained special solutions can be identified with open boundary problems of the associ- 
ated Toda-type systems, and have interesting application to the construction of periodic 

I— —I solutions. 

(N 1 Introduction 
o 

^ The purpose of this paper is to systematise a class of reductions applicable to discrete KdV-type 
QQ equations and introduce a technique to obtain exact expressions for the solution of the reduced 
O systems. The discrete KdV-type equations have a long history going back to [ll|2],[3], cf. and 
^ more recently classification results were obtained by Adler Bobenko and Suris (ABS) O E]. Of 
these systems the type-Q subclass are roughly speaking the least degenerate, they occupy the 
top level of the hierarchy of such equations. The basic reduction we consider is a finite degree-of- 
freedom system which governs the solution on a region of the lattice between imposed singular 
boundaries. Important point for hyperbolic equations, like the discrete KdV-type equations (for 
which there exists a natural Cauchy-problem even when the domain is generalised from a regular 
lattice to the quad-graph [Zl), is that imposed boundary values would normally over-determine 
a solution. It is the singular nature of the boundaries that resolves the expected inconsistency 
and legitimises such reductions. 

This kind of reduction applied to a discrete equation of KdV-type was proposed by Field, 
Nijhoff and Capel in |8] where the lattice potential KdV equation (or HI in ABS classification) 
was considered. The authors linearise the dynamics of the reduced system and give a finite 
iterative procedure to construct its explicit solution (even in the quantum case which is not 
considered here). Similar boundary problems have more commonly been studied for integrable 
equations outside the KdV class. In particular the discrete Hirota equation was considered 



> 

X 



1 



in in] with prescribed singular solution outside of a bounded domain, a situation which there 
was motivated physically. As remarked in [8], probably the nearest related examples come 
from the theory of Toda-type lattice systems [101 E]) where the 'open' (or 'finite non-periodic') 
boundary conditions can be identified with what we call singular boundaries her^ In fact the 
relationship between KdV-type and Toda-type systems, which is well known since Flaschka [12], 
is remarkably nothing but restriction to a sub-lattice for the fully discrete systems [131 [H]. The 
open Toda lattice has a long history going back to Moser [15] and has been extensively studied 
and generalised (see for instance [T6|, [TTl [TBI [191 EOl El]), so is the better-known setting for the 
present work. 

One of our main motives is to understand better the ABS equations as prototypical equations 
whose natural integrability property is the multidimensional consistency. In particular we are 
interested in the primary such model due to Adler |[22j known as Q4 since [5], see [231 E]- 
Remarkably Q4 can be characterised (up to some non-degeneracy) by the symmetry of its defining 
polynomial [23]. On a local level there is an intimate relationship between the singularities and 
the integrability of this equation, in particular analysis of the singular solutions provides an 
explanation of the spectral curve and the natural parameterisation of the equation in terms of 
points on that curve p]. The results reported here stem from the global singularity analysis |25] . 
in particular the solutions obtained can be viewed as soliton-type solutions built on a singular 
background. On the other hand, they can also be viewed as solutions of the open-boundary 
problems for the Toda-type systems connected to Q4, namely discrete analogues of the elliptic 
Toda lattice [2S1EZ] and the elliptic generalisation of the (Relativistic) Ruijsenaars-Toda lattice 

[H ESI EH ED]. 

We proceed as follows. Section |2] recalls the class of equations and their singularity structure, 
and in particular explains the idea of singular-boundary reduction. A particular reduction, 
the singularity-bounded strip, which is chosen for its regularity, is defined more precisely in 
Section |3j The substantial part of the paper. Section |4] through 10, is devoted to obtaining an 
explicit expression for the solution of this reduced system. Section 11 makes connection with the 
reductions obtained by imposed periodicity. The extension of the solution to multidimensions is 



given in Section 12 Conclusions are in Section 13 



2 The equations and singularities in the solutions 

The equations we consider take the form 

Qa,f3{u,U,U,iM) = 0, (2.1) 

were Q is one of the polynomials listed in Table [Tj and a, /3 G C are fixed parameters of the 
equation, u = u{n,m), u = u{n + l,m), u = u{n,m + 1) and u = u{n + l,m -|- 1) are values 
of the dependent variable m as a function of independent variables n, m G Z, therefore equation 



(2.1) connects values of u on the vertices of each quadrilateral of the Z lattice. 

Amongst the similar models listed in [5] we restrict to those in Table [l] because they share a 
common singularity structure [SIES]. Restricting further to the elliptic case, Q4, would also not 



^We do not use the term 'open-boundary reductions' here because it gives the false impression that information 
can pass through the boundary. More or less information passing into the boundary disappears, which is an 
important distinction if we consider the boundary from both sides. 
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Ql^ a{u — u){u — u) — P{u — u){u — u) + aP{a — (5) 
a{u — u){u — u) — (3{u — u){u — u) 

— (3){u + u + u + u — + a(5 — (3"^) 
5 sinh(a)('U'u + uu) — sinh(/3)(-u-u + uu) 

— sinh(a — /3){uu + uu + 6'^ sinh(a;) sinh(/3)) 
Vk sn{a) {uu + uu) — Vk sn(/3) {uu + uu) 

— Vksn{a — /3)[uu + uu — A;sn(Q!) sn(/3)(l + uuuu)] 



Q2 
Q3 
Q4 



Table 1: Canonical forms for the type-Q polynomials [6], the coefficients depend on complex 
parameters a and /3, G C \ { — 1,0, 1} is the modulus of the Jacobi sn function appearing in 
Q4 and for Q3 5 G {0, 1}. 



really lose generality; all similar models can be obtained from it by limiting procedures [HI EI] . 
Nevertheless we find it useful to consider the other systems listed in the table in order to see 
explicit solutions for the full progression through the rational, trigonometric and elliptic models. 
We define singularity following [U]: 



Definition 1. A solution of (2.1) is singular with respect to some vertex of a quadrilateral when 
the equation on that quadrilateral is satisfied independently of the value of the solution on that 
vertex. 

This means that singularities can be characterised by explicit conditions connecting the remain- 
ing three vertices. Solving these conditions leads to another characterisation of singularities 
which is more particular to the equation being considered. For the models listed in Table [T] a 
prominent role is played by an associated function /: 



(2.2) 



Ql' : 


/(O 


= e, 


Q2 : 


fiO 


= e, 


Q3^: 


m) 




Q4: 


no 


= V^sn(^). 



Definition 2. A solution of (2.1) is singular on an edge of the lattice oriented in the n direction 
when on the vertices of that edge it takes the values f{Q and /(C±a) for some choice of sign and 
value of (. Similarly it is singular on an edge oriented in the m direction when on the vertices 
of that edge it takes the values f{() and f{( ± (3) for some choice of sign and value of (. 

The connection between these two definitions is that: 



Lemma 1 (ABS [6J). A solution of (2.1) is singular on some edge of a quadrilateral if and only 



if it is singular with respect to one of the two vertices of the quadrilateral not on that edge. 



The problem of globally describing singularities in a solution of (2.1 ), in the sense of Definition 
is reduced by this lemma to listing the collection of all singular edges in the sense of Definition 
It is also worth remarking that Lemma [l] also has a wider significance: it clarifies the connection 
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(a) 



Figure 1: Admissible singularity configurations on a single quadrilateral in part (a), and attached 
arrows in part (b). The solid lines indicate singular edges. The configurations are given up to 
the symmetries of the square and also up to an overall reversal of arrows in part (b). 

between the singularities and the lattice parameters, i.e. the parameters a and /3 which appear in 
the polynomials of Table [TJ Since the early history these parameters were retained when writing 
such equations in order to encode their integrability, see [SI S] , or more recently for Q4 when 
these parameters enter through the elliptic functions see [221 EH] . 

We call the collection of singular edges the singularity configuration, whilst the set of quadri- 
laterals with singular edges will be called the singular region. They are constrained (generically) 
by some simple geometric criteria: 



Theorem 1 ([25]). Propose a solution of (2.1) with some collection of singular edges which 



determine some singular region. Sufficient conditions for existence of this solution are: 

1. Each quadrilateral within the singular region is one of those in Figure^a) . 

2. It is possible to attach arrows to all singular edges in such a way that (i) each quadrilateral 
within the singular region is one of those shown in Figure^b) and (ii) around any closed 
path of singular edges the arrows sum to the null vector. 

3. The solution on the non-singular region is not overdetermined. 

The conditions 1 and 2, which can be tested geometrically, are sufficient to explicitly construct 
the solution on the vertices connected by the singularity configuration (an example will be given 
in the next section, see also [25]). It is rare that configurations satisfying conditions 1 and 2 
would fail on condition 3, but pathological examples do exist. Also, in non-generic circumstances 
singularities that violate condition 2 can exist, so the conditions of the theorem are not always 
necessary for admissibility. 

Consideration of Theorem [T] reveals that a great variety of singularity configurations are 
admissible. In particular it is easy to specify boundaries which partition the lattice, and if a 
partitioned region has remaining a finite number of degrees of freedom we call the system on 
that region a singular-boundary reduction. Henceforth this paper will focus on one example in 
particular which is introduced in the following section. 
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n 



m 



/(C + 7) 



■/(C-7) 



Figure 2: An admissible singularity configuration which isolates a regular strip of the lattice 
between singular boundaries. Initial data can be placed on the indicated vertices (in the example 
depicted there are M = 6 degrees of freedom). The function / appearing in the solution on the 



singular boundaries depends on which equation is considered and is given in (2.2), the function 
C is defined on as C = Co + '^ct + i^P- Co and 7 are arbitrary constants. 

3 The regular singularity-bounded strip 

Consider the singularity configuration illustrated in Figure |2] which isolates a strip of the lattice 
between singular boundaries. It is clear that initial data on the indicated vertices uniquely 
determines the solution on the strip, the number of these vertices, i.e. the number of degrees 
of freedom, will be denoted by M. This system is probably the most regular singular-boundary 
reduction on 1? allowed by Theorem [!} It is natural to ask what the solution on this strip looks 
like, and answering this question is the basic technical problem we set out to solve in this paper. 
It is useful to formulate a definition to formalise this particular singular-boundary problem: 



Definition 3. A solution of (2.1) is called a singular-boundary reduction on a regular strip with 



M G N degrees of freedom if it satisfies the boundary conditions 



„ = /(C + 7), n + mG {-1,0}, 

W(C-7), n + m G {M + l,M + 2}, > 



and is also defined for all values of n + m E {1, . . . , M}. Here f is the function given in (2.2), 



C = Co + na + m(3, Co O'^d 7 are some constants, and a, (3 are the parameters appearing in the 



equation (2.1) 
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Figure 3: Admissible singularity configurations linked by consecutive Backlund transformations. 
Odd-M case, N e {0,1,2}. 

4 Outline of the solution method 



A central role will be played by the Backlund transformation of (2.1), which takes the form 



(4.1) 



where Q was the polynomial appearing in the equation (2.1) itself. This system determines a 



new solution of (2.1) from an old solution u, the free parameter A* is the Backlund parameter. 
Such Backlund transformations were described in [321 [13] , the lattice and Backlund directions 
are distinguished only by the difference of a parameter. In fact the equation, its Backlund 
transformations, and the superposition principle for pairs of commuting transformations, are 
naturally considered simultaneously as a single system in multidimensions. The compatibility of 
this system is the well-known multidimensional consistency property. 

The key to the solution method is the Backlund chain illustrated in Figure |3} The solutions 
contain singularity-bounded strips of changing width|^ This change in the singularity configu- 
ration is consistent with the Backlund transformation, or in other words, these configurations 
are consecutive two-dimensional cross-sections through an admissible singularity configuration 
in three dimensions. This can be confirmed by applying the conditions of Theorem [T] in Section 
|2} To apply this theorem in three dimensions one simply extends the Definition [2] and calls 
an edge singular in the direction defined by Backlund transformation (4.1) if values on its end 
vertices satisfy u = f{(), = f{( ± X^) for some choice of sign and value of (. 



^The establishment of this sequence is reminiscent of the 'peehng away' procedure applied to the similar 
problem in [5]. 
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Figure 4: Admissible singularity configurations linked by consecutive Backlund transformations. 
Even-M case, N e {0, 1,2}. 

The imposed singularity structure and restricted domain of the solutions are non-generic 
features which will introduce some subtleties into the problem of integrating the Backlund chain. 
The solution in Figure |3]^a) can be identified as the seed solution, but connects only to singularity- 
bounded strips with an odd number of degrees of freedom. A similar sequence in the even-M 
case, illustrated in Figure |4| connects to a different seed solution, therefore with this approach 
it will often be necessary to distinguish between the odd-M and even-M cases. 

From the implied sequences in Figures [3] and |4] the length of the chain, which we denote by 
A^, can be related to the number of degrees of freedom on the strip, M: 

M = 2N-1, Modd, , , 

4.2 

M = 2N -2, M even. ^ ^ 

Generically, application of A^ Backlund transformations introduces 2A^ degrees of freedom from 



the choice of integration constants and Backlund parameters. The apparent mismatch in (4.2) 
is explained because of the constraint Ai + . . . + Aa? = 7 connecting the Backlund parameters 
appearing in the chain (see Figures [s] and |4]) with the fixed parameter 7 appearing in the reduction 
(Figure [2]). In the even-M case the further mismatch is accounted for by the reducibility of the 
Backlund equations connecting the first two solutions of the chain in Figures |4]^a) and|4](b), both 
strips have no degrees of freedom and there is no integration constant involved in the Backlund 
transformation connecting them. 

We remark that the role taken here by the Backlund parameters is revealing, it indicates an 
immediate upshot of this construction: the Backlund scheme reduces what would be an M*'*- 
order integration of the system on the strip, to just A^ = (M + l)/2 first-order integrations if M 
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is odd, and N — 1 = M/2 first-order integrations if M is even (A^ — 1 instead of because as 
noted above the first Backlund transformation trivialises in the even-M case). 



5 Integrating the Backlund scheme 

The problem of integrating a Backlund chain of length for this class of equations was reduced 
in [33] to solving a set of N independent first-order, homogeneous, linear, but multidimensional 
equations. The method relies on the existence of 2^ particular solutions of the Backlund chain 
which arise naturally by assuming that the seed solution itself is multidimensional. This section 
is devoted to recalling that result which will play a central role in construction of the singularity- 
bounded strip solution described in Definition [3j 

Independent variables in + 2 dimensions will be denoted n,m,li . . . In & the correspond- 
ing shift operators and lattice parameters are ~, " , Ti . . . Ttv and a, /3, Ai . . . Xn respectively. The 



multidimensional system generalising (2.1) is simply 



Qa,X. 



U,U,JiU,JiU) 



0, Q/3,A,(M,M,TiM,TiM) =0, i G J, 



(5.1) 



Qx„x, (u, JiU, JjU, JiJju) = 0, i, j G J, 
for some J C {1 . . . A^}, which is a system for a function u defined on Z^"*"^, but we have allowed 



for J to perhaps be only some subset of indices. Likewise, the Backlund transformation of (5.1) 



is obtained by extending (4.1) 



QA„A.(«,TjM,M=,,TiMj 



: 0, i e J, 



(5.2) 



which takes a function u defined on Z^+^ satisfying (|5.1|) and defines a new solution u^,. 



The following theorem applies to the integration of a Backlund chain when the seed solution 



ti of (2.1) on 7? has some known extension satisfying (5.1) on Z^"*"^. 



Theorem 2 ([S3]). Let u : — ^ C U {oo} he a solution of (^j]) in the case J = {I...N}. 
Suppose for each i G {1 . . . A^} that function rji : Z^+^ — > C U {oo} is such that defined on 
expression 

' ^'""-^^ (5.3) 



satisfies (5.2) with J = {1 . . . N} \ {i}. And define Mi...Ar by expression 



[l-0iT?]...[l-0^T^]t; 



.Till' 



ie{i...Af} 



[n 



N 

je{l...N}\{i} 



TJ']v^ (5.4) 



[1 - (piTj] ■ ■ • [1 - (/^TV I ]vJ 

(here v and 0i . . . o,re just shifted versions of u and rji. . .rj^)- Then jv|ii=.. -i^=o is a 



solution of (2.1) related to , =/j^=o by the composition of N Backlund transformations (4-1) 
with Backlund parameters Xi . . . X^. 
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This theorem decouples what would be an iterative procedure to integrate the Backlund chain 
into integration of N separate systems, one for each of the auxiliary variables rji . . .tjn. Further- 
more it was shown in [33] that substitution of (5.3) into (5.2) with J = {1 . . . N} \ {i} yields a 
homogeneous linear system for each rji. 

The additional property of rji...rj^, that (Tjr7j)(Tj^?7j) = (Tj?7j)(T~^?7j) Vi,j G {l...A^} 
which was also proven in |33], implies existence of a potential for these variables: 



Lemma 2. Suppose u is a solution of (5.1) with auxiliary functions t]i 
Theorem Then there exists a function r defined on Z^+^ such that 



iE{l...N}. 



rjN- as described in 



(5.5) 



Written in terms of this potential (5.4) acquires the simpler form 



UI...N 



[Tr' - Ti] • ■ ■ [TJ,' - TN]ru 
[Tr^ - Ti] ■ ■ ■ [J-, 



TN]r 



(5.6) 



But in fact our motive for modifying Theorem [2j by introducing r and (5.6) is not just cosmetic, 
it will allow to incorporate singularity structure of mi...7v which is not possible through rji . . .rjj^ 



and (5.4). 



Our objective is to apply the above (modified) multidimensional integration procedure to the 
Backlund chain with singularity structure described in Section |4j Thus it necessary to consider 
not just the first and last, but also all intermediary solutions along the chain. Solutions resulting 
from application of transformations with parameters {Aj : z G /} for any subset of indices 
/ C {1 . . . A^} are implicit in (5.6) by replacing the indices {1 . . . A^} with this subselj^ 



Ui 



[n..,[Tr^-T, 



r 



I (Z{l...N}. 



(5.7) 



This expression for the intermediary solutions will play the central role in the following two 
sections. First it allows conversion of the imposed singularity structure into constraints on 
r. Second it transforms the Backlund equations along the chain into (multiplicatively) linear 
equations for r. 



6 The singular boundaries in multidimensions 

This section is devoted to establishing preliminary constraints on the multidimensional potential 



r such that the expression for the solutions along the Backlund chain, ui in (5.7), satisfies the 
desired singular boundary conditions. 

The boundary conditions are the realisation of sequences implied by Figures |3] and |4] in 
multidimensions. They are central to the solution procedure so we introduce some terminology 
for them: 



•^We adopt the notational convention that U{} = u, u^i^ = u.i, u^i j^ = u.ij etc., i.e. we omit the enclosing 
braces and separating commas on the set of indices of the sohitions along the Backlund chain. 
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Definition 4. For fixed integers M and N satisfying (4-2) and set of constants Ai . . . A^r G C 
we call the functions Uj, / C {1 . . . A^}, an epitaxial sequence if they are defined on the subset of 



n + m + h + ... + In e {N 
and satisfy the conditions 



Ul 



fiC + E.e/ A 



n + m + li + . . . + In 
n + m + li + . . . + In 



G 



{N- 
{M 



M- |/| +2} 



2}, 



(6.1) 



on the boundary of this region. Here 
function 

C = C,Q + na + m(5 + ZiAi + . . . + In^n 
defined on Z^+^, C,q is some constant. 



N +\I\ + l,M - N +\I\ 
denotes the number of elements in I, and ( is the 

(6.2) 



The epitaxial sequence generalises to multidimensions and to generic / C {l...A^} the 
singularity configurations in two dimensions for I = {}, / = {1} and / = {1,2} which were 
illustrated in Figures |3] and |4j Conditions (6.1) reduce to (3.1) in the case / = {l...A^} and 
/i = . . . = /at = (corresponding to the last solution of the chain restricted to Z^). And in the 
case / = {} (corresponding to the first solution of the chain) they become 

/(C), n + m + h + + e {N -1,...,M - N + 2}, 



u 



(6.3) 



completely defining the seed solution, m, which participates in the expression for uj (5.7). Note 
that the multidimensional strip on which the seed solution is defined by (6.3) has width 3 if M 
is odd and width 2 if M is even (cf. (4.2)), and we do not define it outside this strip. It is useful 
to introduce terminology for this region: 



Definition 5. For fixed integers M and N satisfying (4-2) the subset ofTl^^"^ where the first 
function of the epitaxial sequence is defined, i. e., where n + m + li + . . . + Ij^ G {A^ — 1, . . . , M — 
A^ + 2}, is called the seed domain. 

There is now an elegant thing that happens, captured in the following lemma: 



Lemma 3. Expression (5.7) defines an epitaxial sequence corresponding to integers M and N 
satisfying (4-S) and constants Ai . . . Xn if and only if u is defined on the seed domain by (6.S) 
in terms of ( (6.2), and the associated function r satisfies 

r = 0, n + m + /i + ... + /iv G {-A^- 1,...,M + A^ + 2} \ {A^- 1,...,M- A^ + 2}, (6.4) 

and is also defined throughout the seed domain. 

Essentially r is required to vanish on the region that the seed solution u is un-defined, regions 
of the lattice further beyond where r vanishes in (6.4) do not participate in the problem. 

Lemma |3] may be proven as follows. The first statement of the lemma is obvious, it cor- 
responds to the case / = {}. That (6.4) is sufficient for ui defined by (5.7) to satisfy (6.1 



can be verified directly by calculation. We give a simple recursive procedure to show that (6.4) 



is necessary for ui defined by (5.7) to satisfy (6.1). Assume that (6.1) holds for all subsets 
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of indices / C {l...A^}. For some / C {1 . . . A^} suppose that r vanishes on the domain 
{N - 2\I\ -1,...,M-N + 2\I\ +2}\{N -1, . . . , M - N + 2} (importantly this is trivially true 
when / = {}). Then by calculation show that if ( 6.1[ ) is to hold when /—)■/□ {i} {i ^ I), then r 
must also vanish on the region {A^ - 2|/| - 3, iV - 2/ - 2} U {M - + 2| J| +3, M- A^ + 2|/| +4}. 
Apply this fact iteratively to the sequence / = {}, / = {1}, / = {1,2}, . . . to conclude that ( |6.4 ) 
must hold. 

Thus in Definition |4] we have formulated the desired multidimensional singularity structure 
of solutions along the Backlund chain, which we have called an epitaxial sequence. Lemma [3] 
shows that this required singularity structure translates through (5.7) into the remarkably simple 
constraint that r should vanish wherever the seed solution is un-defined. In the following section 
equations complementing (6.4) will be established that determine r also where the seed solution 
is defined. 



7 Equations for r on the seed domain 

This section is devoted to establishing equations for r resulting from the requirement that uj in 



(5.7) satisfy the equations along the Backlund chain. As the theorem in Section 5 implies, the 



equations linking the first two solutions of the chain would usually be sufficient to determine r 
through auxiliary functions rji . . .7]^. However the situation here is non-generic, specifically the 



already established constraint (6.4) means that within the required domain these equations are 



in fact already satisfied. Instead we need to go further along the chain, exactly how far depends 
on whether M is odd or even. The result of this section is summarised in Lemma H] at the end. 

In the odd-M case primitive equations for r within the seed domain, i.e. the region in 
where n + m + h + . . . + In G {A^ — 1, A^, A^ + 1}, can be obtained by requiring u,, for generic 
i E {1 . . . A^}, satisfy with J 
u Ui for generic i G {1 . . . A^}, 
equation in the m direction: 



{1 . . . A^} \ {i}. Consider the first equation of (5.1) with 
but, for the sake of symmetry later on, backward-shift the 



Ui] 



0. 



(7.1) 



For brevity we consider only the first equation, the other equations in (5.1) can be obtained 



from this one by permuting the lattice directions. Given u in (6.3) and the established vanishing 
of T in (6.4), the terms appearing in (7.1) obtained from (5.7) evaluated when / = {%} and 
n + m + li + ... + l]<q = N are 



T-^f)/(C-A,)-(T,f)/(C + A. 



Ui 



T-V)/(C-A.)-(T,r)/(C + A.) 



(TrV) - (T,;r) 



Ui = /(C - Aj)- 



(7.2) 



Substitution of these expressions into ( |7.1 ) results in the following equation for r: 



(Ta(TrV) _ QaAfiC + A.), /(C - A.), /(C + A,), /(C - A.)) 
QaAf iC - Ai), fCC - A.), /(C + A.), /(C + A.)) ' 
n + m + li + . . . + In = N, M odd. 



Trf)(T.r) 



(7.3) 
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This can be verified using the fact that the polynomial Q is degree-one in each variable and that 
(7.1) has particular solutions Ui = Tf^u (the same method used to obtain equations for rji in 
[33J). The equation (7.3) connects values of r on the domain n+m+li + . . . + In G {-^ — 1, A^+1}. 
By construction, the equation ( |7.3 ), plus similar equations obtained from this by permuting the 
indices and lattice directions (i.e., resulting from other equations in (5.1)) are necessary and 
sufficient conditions on r for the expression in (5.7) to satisfy the Backlund chain of length one. 
The result of [33] can then be applied to conclude that expression (5.7) also satisfies the chain 
of arbitrary length. 

In the even-M case the primitive equations for r on the seed domain, which is the strip 
n + m + li + . . . + Ifyf G {A^ — 1, A^}, can be obtained from the Backlund equations relating Ui 
and Uij for generic i,j G {1 . . . A^}. Specifically (5.2) in the case J = {1 . . . A^} \ {i,j}, u — )■ Ui, 
and A* — )■ A^, so for instance the first equation of (|5.2[) becomes 



Qa. 



0. 



(7.4) 



Given u in (6.3) and the established vanishing of r in (6.4), now in the even-M case, the terms 
in (7.4) obtained from (5.7) evaluated when n + m + li + ... + li\i = N — 1 are 



Ui 



/(C - A. 



Ui 



Ui 



(TriT,r)/(C - A. + A,) + (T.TtV)/(C + A. - A,) 

{J;'J,t) + (T,TtV) 
(T-^T,f)/(C - A, + A,) + (T.TTif )/(C + A. - A,) 



(7.5) 
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-T,r) 



(T.TT 



The substitution of these expressions into (7.4) yields the equation 



(T.TT 



1^ 



f)(T 



-1 



T,r) 



Qa,x.ifiC + A,), fiC - A,), /(T.C - A,), f(J-'C + A,)) 

(TriT,f)(T,T-V) Q^,xXf{C + A,), /(C - A,), /(T'^C + A,), /(T,C - A,)) ' 

n + m + li + ... + lN = N-l, M even. 



(7.6) 



which connects values of r within the desired domain n + m + h + . . . + In G {A^ — 1, A^}- In 
this case the equation (7.6), plus the equations obtained from this by permutation of the indices 
and lattice directions, are necessary and sufficient conditions on r for the expression in (5.7) to 
satisfy the Backlund chain of length two. This chain is long enough for the solution to emerge 
beyond the singular situation (i.e. a sufficiently large value of A^ for M to be greater than zero), 
and therefore the expression (5.7) also satisfies the chain of arbitrary length [33j . 

The multiphcatively linear equations (7.3) and (7.6), plus similar equations obtained by 
permuting indices and lattice directions, are the desired constraints determining r on the required 
domain in the odd-M and even-M cases respectively. The expressions on the right-hand-side of 
these equations can be evaluated using the explicit forms of the polynomials Qa,i3 in Table [l] and 
the corresponding functions / in (2.2). Details of this calculation will be given in the following 
section, for now we give the results. In all cases they reduce to the following autonomous 
equations for an associated variable vr: 



(T,f)(TrV) g{P)gia-Xi 



(Tr'^)(T,7r 



T,7r) 



g{a)g{l3 - Ai) 
g{Xj)g{a - Aj 



(TriT,-^)(T,TTV) (7(A,)^7(a - A,; 



n + m + li + + = N, M odd, (7.7) 
n + m + li + ... + lN = N-l, Meven, (7.8) 
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which involve a new function g dependent on the equation in question, 



Q1\Q2 : 


9(0 


= e, 


Q3' : 


9(0 


= sinh(^), 


Q4: 


9(0 


= H(Oe(0, 



(7.9) 

and where the variable n is related to r differently in the case of QA than for the others: 

n + m + h + ... + In e {N -1,...,M - N + 2}. (7.10) 



Q1\Q2,Q3^ : 
Q4 : 



r = TT, 
r = 0(C)7r, 



Thus the problem of constructing r on the seed domain is reduced to obtaining the solution of 



autonomous systems for vr formed from (7.7) or (7.8) by permuting indices and lattice directions 



however these systems for vr can be further improved. 

A-priori vr need only be defined when n + m + /i + ... + Zjv £ {N — l,...,M — N + 2} (i.e. on 
the seed domain), however relaxation of the domain of vr to the whole of Z^"*"^ turns out to be 
quite natural. The reason is, this permits the systems for n, in both odd-M and even-M cases, 
to be replaced by the one simpler system 



T,^)(T7i7r) _ g{a-X,) {T^iT){Tr^7r) _ g{/3 - X,) 

9{m^^y 



(Trifr)(T,7r) g{a)g{X.y {J^'7^){J,'n 
(T,T,vr)(T-M_^?(A,-A.) 



zG{l...iV}, 



(7.11) 



(TriT,vr)(T,7r) g{X,)g{X,y 



It can be directly verified that when vr is restricted to the domain n + m + /i + . . . + / 



N 



{N — 1, + 1} the system (7.11 ) is equivalent to the system formed from (7.7) by permuting the 



indices and lattice directions. Likewise, the system formed from (7.8) by permuting the indices 
and lattice directions is equivalent to (7.11) restricted to the domain n + m + li + . . . + In G 
{A^ — 1, A^}. Thus, regardless of the value of M, vr is without loss of generality governed by 
7.11) on Z^+^. Note that consistency of (7.11) relies only on the basic feature of functions in 
T^^that g{-0 = -g{0. 

The following lemma delivers the result of this section. 



Lemma 4. Fix integers M and N satisfying (4 -2} and constants Ai . . . Aat G C. Suppose u and t 
satisfy the constraints of Lemma^ and Uj, / C {1 . . . A^} is the corresponding epitaxial sequence 
(Definition^ obtained from u and t by (5.1). Furthermore suppose that on the seed domain r 
is of the form (7.10) where vr defined on Tj"^ satisfies ( 7.11). Then for generic / C {1 . . . A^} 
and i G {1 . . . A^} \ / the substitutions u ^ uj, — )■ "U/uj,} and A* — i- Aj satisfy the Bdcklund 
equations (5.2) with J = {1 . . . A^} \ / \ {i} throughout the region on which uj is defined (cf. 
Definition |^pL 



Loosely speaking, we have obtained sufficient conditions on r for the epitaxial sequence in (5.7) 



to also be a Backlund chain. It means the integration part of the problem has been captured in 



the multiplicativly linear, autonomous, multidimensional system for vr (7.11). The calculations 
within this section establish the validity of Lemma |4] except for the single ommission that was 
mentioned, which we turn to now. 
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8 Simplification of the equations for r: elliptic case 



The simplification of equations (7.3) and ( |7.6[ ) governing r resulting in equations (7.7) and (7.8) 
for tlie associated variable vr can be confirmed directly by calculation. However, in the elliptic 
case some ways of performing the calculation are much easier than others, therefore this section 
is devoted to giving a set of intermediary steps which are easy to verify. Definitions of the Jacobi 
elliptic and theta functions involved in the calculation can be found in Chapter 5 of [3l] or 
Chapter X of |aSj. 

In this section Qa^/s will denote the Q4 polynomial appearing in Table [T| and the associated 
function / (2.2) is taken to be 

/(0:=v/fcsn(0 = |||. (8.1) 

Central are the following basic identities involving the Q4 polynomial on a quadrilateral with a 
singular edge 



Qa,tifiO,fi^ + a),y,z) 



-- /(a)t(e, b,a-b) [/(e + b)-y] [/(^ + a - b) - z], 
/(a)t(e, -b, b - a)[/(e -b)- y][f{^ -a + b)-z], 



where we have introduced the symmetric function 



t(a, 6, c) 



l + f{a)fib)fic)fia + b + c) 



Q{0)Q{a + b)e{b + c)e{c + a) 

e(a)e(6)e(c)e(a + 6 + c) 



(8.2) 
(8.3) 



Identities (8.2) and (8.3) can be verified directly from the Q4 polynomial as given in Table [l| 
using the addition formula for the sn function and the first expression for t in (8.4). Alternatively 
these identities can be verified using t in its second guise, whilst starting from the three-leg-form 
of Q4 [SIEB]: 



■X 



Qa,bif iO, 2 _ ^) Q(^) 0(5) Q(^) + + K) 

([/(e + a) - x][/(e - 6) - y]\f{^ _ a + 6) - e(e + «) -b)Q{^-a + 6)- 

[/(e -a)-x] [f{i + h)-y][f{i + a-b)-z\Q{i-a)Q{^ + b)Q{^ + a-h)), (8.5) 

where K is the standard notation for the quarter-period of sn function satisfying sn(i^') = 
1, sn'(i^) = 0. This requires the following basic identity expressing the difference of elliptic 
functions ratio of products of theta functions: 



/(e + a)-/(e 



H(a)e(a)H(^ + i^)e(e + ^) 



Q{i + a)Q{i-a)YL{K)Q{Ky 
The basic symmetries of the canonical type-Q polynomials |14j : 

Qa,b{w,x,y,z) = Qa,a-b{x,w,y,z) = Qh-a,b{y,x,w,z) = -Qb,aiz,x,y,w), 



(8.6) 



(8.7) 



will be used to bring the expressions on the right-hand-side of equations (7.3) and (7.6) to a 
form in which the identities (8.2) and (8.3) can be directly applied, the rest of the calculation is 
then quite straightforward. 
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Starting from expression (7.3) there follows the string of equalities 



(T,f)(TrV) QaAf iC + A.), fCC - A.), /(C + A.), fCC - A.)) 



(T.-'f )(T.r) Qo^Af iC - A.), fCC - A.), /(C + A.), fCC + A.)) ' 
_ QUfiC + A.), /(C + A.), /(C - A.), /(C - A.)) 
QaAf iC - A.), /(C - A.), /(C + A.), /("c + A,)) ' 
_ fimC + -c^, « - mi( + /(C - A. + a)] 

/(a)t(C - A„ /3, « - /3) [/(C + A, - /3) - /(C - A, + /?)] ' 
fimC + A^, -a, a - /3) e(C - A. + /3) e(C + A. - /3) H(a - 



(8.8) 



/(a)t(C - A„ /3, a - /3) e(C + A^ - a) e(C - Xi + a) H(/3 - A^) e(/3 
e(C + A» + g - /3) e(C - A^) H(/3) e(/3) H(a - A^) Qja - A^) 
e(C - A, + a - /3) e(C + Ai) H(a) e(a) H(/3 - A^) e(/3 - A^) ' 



A.) 



Here the symmetry Qa,f,{w , x , y , z) = — Qfe,a(t(7, x, z) has been applied to the numerator of 
the first expression on the right-hand-side in order to obtain the second expression. The third 



expression follows from the second by substituting for shifts on ( using (6.2) and then applying 



identities (8.3) and (8.2) to the numerator and denominator respectively. The fourth expression 



results from the third by applying identity (8.6) to the factors appearing in both numerator 



and denominator. The fifth expression emerges after substituting for / and t from (8.1) and 
(8.4). This final expression motivates the substitution r = <d{()n (7.10) and identification 



g{0 = 0(0 H(0 ('''•9) resulting in equation (7.7) governing the new variable vr. 



A similar string of equalities can be obtained starting now from expression (7.6): 



(T,T7'f 



)ar^T,r) 



Qa,xXf iC + A,), fiC - A,), /(T.C - A,), fjJj'C + A,)) 

(Tr%.f)(T,T-V) Qa,xAf iC + A,), /(C - A,), /(T'^C + A,), /(T,C - A,)) ' 
_ Qa,^..a,(/(T.C - A,), /(C - A,), /(C + A,), /(T-^C + A,)) 
Qx.Afi'^i'C + A,), /(C + A,), /(T,C - A,), /(C - A,)) ' 
/(A, - a)ti( - A, + \, -K, a) [f{( - A,) - /(C + A,)] 
/(A,)t(C + A,- - A„ a, K - a) [/(C - A, + a) - /(C + A, - «)] ' 
_ f{Xi - a)t{C - Xj + A„ -Ai, a) e(C + Xj - a) e(C - Xj + a) H(Aj) Q{Xj) 



(8.9) 



/(A,)t(C + Xj - A„ a, Xi - a) e(C + A,) e(C - Xj) H(A, - a) Q{X, - a) ' 
^ e(C - Xj + Xi + a) 9(C + Aj - A^) H(A,) Q{Xj) H(A^ - a) e(A, - a) 
e(C + Aj ~Xi + a) 0(C - A, + Xi) H(Ai) Q{Xi) H(A, - a) Q{Xj - a) ' 

Here the second expression is obtained by application of the symmetries Qa,b{w,x,y, z) = 
Qb-a,b{y,x,w, z) and Qa,b{w , X , y , z) = —QbAu^'^^^T^) to the numerator and denominator re- 
spectively of the first expression. The third expression is obtained from the second by substi- 



tuting for shifts on C, using (6.2) and then applying identities (8.3) and (8.2) to the numerator 



and denominator. The fourth expression is again obtained from the third using (8.6) and the 



fifth expression emerges from the fourth after substituting for / and t. The final expression here 



motivates the same substitution as before, written in (7.10), and results in (7.8). 
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9 Plane-wave factors 



It was described in Section[5]how in the generic case the last solution of a Backlund chain, ui,,,n, 
can be expressed in terms of the seed solution, u, and either r ( |5.6 ) or the variables 7]i. . .T]]y 
(5.4). The solution procedure followed in Sections |6] and [T] neatly handles the imposed singularity 
structure by focusing on the potential r, but now we ask about rji . . .rj^. 

Consideration of the constraint on r (6.4) which resulted from the imposed singularity struc- 
ture shows that associated variables rji . . .1]^ remain un-defined by relation (5.5) on some partic- 
ipating parts of the lattice. However the relaxed domain for the associated variable vr introduced 
at the end of Section allows introduction of variables defined on Z^"*"^ 



as 



T,;7r 



Qi 



TJ 



TT 



ie{l...N}. 



(9.10) 



Occurences of vr in the solution Mi...Ar can be expressed purely in terms of variables gi . . . qn ii^ 
almost exactly the same way as leads to its expression in rji . . .rji^ for the non-singular case. 
Furthermore the system (7.11) governing tt can also be written just in terms of these variables, 
whence it takes the form: 



Qi 



g{a)g{Xi 

g{^j)g{^i) 



Qi ) Qi 



-Qi, 



^ gi(3 - A, ) 

gmgi>^i 
t,je{i...N}, 



zE{l...N}, 

3 7^ 



(9.11) 



One advantage of these variables is that (9.11) can be integrated by inspection: 

' g{a - Xi 



Qi 



gW - A.) 
g{a)g{Xi)J \g{l3)g{Xi) 



n 

je{i...N}\{i} 



^(Aj - Xj 

g{^j)g{^i 



(9.12) 



for each i G {1 . . . A^}, where q is an arbitrary function of the single variable Zj. Establishing the 
corresponding solution of (7.11) requires more effort, it looks as follows 

'g{a - Xi) 



-= n 

ie{i...Af} 



g{(^)g{K 



gW - A,) 
gimK) 



n 

je{i...N},j<i 



g{.^j - K 
g{^j)g{^t 



Y[ (-i)'^('' 

je{i...N},j>i 



-1) 



k/2 



(9.13) 

where for each i G {1 . . . A^} c* is an arbitrary function of U (which is related to Ci above). 

This completes the integration part of the problem, it remains to write down the solution by 
substituting the explicit expressions for u and r in (5.6). 



10 Explicit form of the solution 



The constructed multidimensional potential r (cf. Lemmas and and expression (9.13)) yields 



through (|5.6|) a representation of the regular singularity-bounded strip solution, specifically the 



restriction to of 



UI...N 



Ti] ■ ■ ■ [T^^ 



TA^]r/(C) 



[Tr^-Ti]---[T^^-T 



NT 



(10.1) 
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where / was given in (2.2) and ( in (6.2). To evaluate this expression in terms of functions 
defined only on requires addressing some final combinatorial issues due to the vanishing of r 
(pl). 



We require to expand the operator present in the numerator and denominator of (10.1), 
consider first the denominator: 



IC{l...N}^ 



i&{l...N}\I 



10.2) 



Due to (6.4) terms in this expansion which evaluate r outside the seed domain (cf. Definition 



5j) will be zero, more precisely, for the expression 



+m+ii+...+i]v=M 



:io.3) 



to be non-zero requires /i + |/| — (A^ — |/|) G {N — 1, . . . , M — N + 2}. Explicitly as a condition 
on |/| this reads 2|/| G {2N — 1— /x, ...,M + 2 — //}. Analysis of this condition by separation 
of the cases when M and /i are odd or even leads to the following: 



Lemma 5. Suppose for some integers M and N satisfying (4-S) that t satisfies the constraint 
{6.4) ■ Then for expression (10.3) with fi G { — 1, . . . , M + 2} to be non-zero it is necessary that 



\I\ G Gf^ where 



{(M-/i)/2,(M-^ + 2)/2}, 
{(M-^ + l)/2}, 
{(M-^ + l)/2}, 
{(M-/. + 2)/2}, 



odd-fi, odd-M, 
even-fi, odd-M, 
odd-fi, even-M, 
even-fi, even-M, 



fie {-1,...,M + 2}. (10.4) 



It is useful to see G^ explicitly for small values of M, choosing M G {0, 1,2,3,4,5} the vector 
[G-i . . . Gm+2]'^ looks as follows: 



{1} 
{1} 
{0} 
{0} 



{1,2} 

{1} 
{0,1} 
{0} 
{-1,0} 



{2} 
{2} 
{1} 
{1} 
{0} 
{0} 



{2,3} 

{2} 
{1,2} 

{1} 
{0,1} 

{0} 
{-1,0} 



{3} 
{3} 
{2} 
{2} 
{1} 
{1} 
{0} 
{0} 



{3,4} 

{3} 
{2,3} 

{2} 
{1,2} 

{1} 
{0,1} 

{0} 
{-1,0} 



(10.5) 



Of course which terms are non-zero in the expansion of the numerator of ( jlO.l ), where the same 
operator acts on some function multiplied by r, is also determined by the condition |/| G G^. 

Evaluation of the non-zero terms in the expansion (10.2) is straightforward given (7.10), it 
comes down to the following identity 



■i)'^'[n...T.][n 



ie{i...N}\i 



vr 



[n 



ie{i...N} 



TT 



Yi n.e/ ( [n,e{i...iV}\«T/]^?' 



IC{1...N}, 

(10.6) 
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where the denominator here is a normahsing factor which may be given to all terms in the 
expansion, and we have introduced important constants Yj for each J C {1 . . . A^} defined as 



Y, -■- 



■1)"' n 



I C{1...N}. 



;io.7) 



The identity (10.6) follows directly from (9.10) and (9.11). Apart from constants Yj the other 
objects appearing in (10.6) are the plane-wave factors Qi . . . Qn, but to express the solution we 
are interested only in their restriction to Z^, it is therefore convenient to introduce functions 
defined on to replace them, we define 



Pi 



Pifi 



9{(^ - 
g{a)g{Xi 



9W - X^) 
9m9{^^) 



ie{l...N}, 



10.8) 



where the constants pi^ . . . pn,o come to replace ci . . .cn appearing in (9.12) and are chosen so 
that 



=«jv=0 



(10.9) 



Expanding ( 10. 1[ ) then using Lemmajsj identity (10.6) and relation (10.9) yields the following 
more explicit form of the solution: 

Theorem 3. To each polynomial in Table [I] associate complex functions of a single complex 
variable x, y and g as follows: 







y{0 


9{0 


Q2 

gs^ 


e 

H(0 


1 
1 

1 

e(0 


sinh(^) 



(10.10) 



Choose M en, set N = {M + l)/2 if M IS odd, and N = {M + 2)/2 if M is even, and fi x some 
set of parameters Ai . . . A^v G C. Define constants Yj for each / C {1 . . . A^} by (10.1), define 
functions pi. . . on I? by (10.8), and define sets G_ 
{(n, m) GZ^:n + mG{— 1,. 



Gm+2 by (10.4). Then u defined on 



u 



M + 2}} by the expression 
Z]/c{i...Af},|/|eG„+„^^ (riie/Z^O ^(C + Yliei 



ie{i...N}\i 



A,;) 



E 



/c{i...Af},|/|eG„+, 



i&{l...N}\I 



(10.11) 



is a solution of (2.1) satisfying Definition^ the function ( appearing in (10.11) coincides with ( 
from the definition, whilst the constant 7 appearing in the definition is given by = Xi + . . . + Xn . 



The obtained expression ( 10.11 ) concludes our construction of the regular singularity-bounded 
strip solution on Z^. Degenerate solutions can be obtained from (10.11 ) in the limit as the differ- 
ence between some of the parameters Ai . . . Ajv approach either zero or a half-period of function /. 
The method of construction allows us to conclude, except for these degenerate cases, that (10.11) 
is the general solution satisfying Definition |3| There remains of course the inverse problem, to 
obtain pi^o • • • Pn,o and Ai . . . Aat in terms of some Cauchy data, say m(0, 1) . . . ^(0, M). 
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Although the multidimensional origin of this solution is more apparent through the repre- 
sentation (10.1), the fully explicit form of the solution (10.11) exposes important features. The 
solution is in essence a ratio of Hirota-type polynomials; for the class of equations being con- 
sidered this is a generally expected feature of solutions obtained through Backlund iteration. 
The most prominent new feature of expression (10.11) is that the Hirota-type polynomial is 
segmented across the width of the strip. Specifically, at some fixed value of n -|- m m is a ratio 
of Hirota-type polynomials which contain only terms homogeneous in pi. . .p^ of homogeneity 
degree given by Gn+m (10.4), (10.5). Despite this segmentation there are clearly identifiable 
aspects of the Hirota-type polynomial, such as constants Y^ijy, which are new and distinct from 
the analagous quantities present in the soliton-type solutions of these equations obtained in 

[STlEHlEniliniES]. 



11 Relation to periodic solutions 

The explicit solutions constructed here generally need not be defined outside of their natural 
domain, which is a diagonal strip of the lattice (cf. Section |3|. However, strips could of course 
be glued together, in principle obtaining a solution that fills the lattice. 

It is interesting to consider the result of also imposing periodicity on such constructed solu- 
tions. The so-called pg-reduction for some fixed non-negative integers p and q imposes that 

u{n + p,m + q) = u{n,m), n,m&'Ij, (H-l) 

and to avoid trivialisation we also assume p and q are not both zero. Such periodic problems 
for discrete KdV-type equations have received a fair amount of attention, the earliest work 
establishing basic techniques based on the integrability is due to Papageorgiou, Nijhoff and 
Capel [H], and since then the problem has generated a substantial literature, see for instance 
[I2l HSl EH SHI SHI 113 EH]- Note that assuming p and q are non- negative here does not lose 
generality due to reflection symmetry of the considered systems. 

The most basic system constructed by imposing both periodicity and singular-boundaries is 
probably to require a solution of (2.1) satisfy Definition [s] plus ( |11.1[ ), whilst choosing 



M + 2 = p + q, 2-f=pa + q/3, (11.2) 

to ensure consistency of the two reductions. Effectively we glue repeatedly the same singularity- 
bounded strip solution in order to fill the lattice. This situation is depicted in Figure |5| 

Such gluing procedure constructs a periodic solution from the singularity-bounded strip so- 



lutions. The additional constraints (11.2) fix 7 and M, leaving all other free parameters of the 
solution un-constrained. The M degrees of freedom on the strip plus the remaining parameter of 
the reduction (q add up to a total of M + 1 free parameters in the periodic solution constructed 
this way. This is just one parameter less than the p + q = M + 2 parameters that should be 



present in the most general solution of the periodic problem (11.1). It poses the question of 
how these singular-boundary reductions fit into the wider theory of the periodic problem for 
this class of equations. Especially whether they have a role to play in obtaining more explicit 



parameterisation of the general solution satisfying (11.1) 
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Figure 5: Imposing the pq periodic reduction simultaneously with the singular boundaries. The 
value of the solution coincides on the two vertices indicated by filled circles, initial data for the 
problem can be placed on the vertices with open circles. 

12 Extension to multidimensions 

Although multidimensional consistency plays the central role in the solution method, the actual 
singular-boundary problem (cf. Definitional) has been considered so- far only on Z^. Solutions 



on Z for arbitrary d eN are also allowed, whence (2.1) is generalised to the system 



Qa,,a,{u,SiU,SjU,SiSjU) =0, z, j G { 1 . . . , (12.1) 

and where we have introduced shift operators Si ... and lattice parameters cti . . . replacing 
~, ' and a, f3, whilst u = u{ni, . . . , Ud) for corresponding independent variables ni . . . replacing 
n and m. This section is devoted to generalisation of the regular singularity-bounded strip, and 
its solution, to this higher dimensional setting. This is motivated by the works [T3| 1^ : 
some related systems of interest are obtained by restriction to a sub-lattice in Z"^ with d > 2. 
We make the higher- dimensional extension explicit here because the solutions are not defined 
isotropically, so this kind of generalisation might not be so immediately obvious. 
The natural extension of Definition [3] is 



Definition 6. A solution of (12.1) is called a singular-boundary reduction on a regular strip 



with M eN degrees of freedom if it satisfies the boundary conditions 

„_/ /(C + 7), n, + ... + ndE{-l,0}, , . 

l/(C-7), ni + ... + n,G{M + l,M + 2}, ^'^'^^ 

and is also defined for all values of ni + . . . + rid E {1, • • • , M}. Here f is the function given in 
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(2.2), C = Co + niai + . . . + naad, Co o,n'd 7 are some constants, and ai . . .ad are the parameters 
appearing in the system (12.1). 

This allows formulation of the following: 

Corollary (to Theorem [s]) . With the same suppositions as Theorem^ except for the definition 
of functions pi . . . pjy ihat are to be replaced by functions defined on 7r, 

Pr.= P.J'-^r^T('-P^T---('-^r^ (12.3) 



^g{ai)g{Xi)J \g{a2)g{Xi)J \g{ad)g{Xi 
the function u defined on {(jt-i, . . . , Ud) & 7/ : ui + . . . + Ud E { — 1, . . . , M + 2}} by the expression 

^ Y,IC{l...N},\I\€G^^I (riigj Pi) ^(C + ZligJ ~ Y.ie{l...N}\I -^») (12 4) 

Z]/C{l...Af},|/|eG^^^ (Hie/ Pi) + -^i ~ J2ie{l...N}\I '^i) 

where jj = ni + . . .+nd, is a solution of (12.1) satisfying Definition^ the function C appearing in 
(12. 4) coincides with the one in the definition, whilst the constant'^ is given by = \i + . . . + \n . 

The generic (finite and with no self-intersecting characteristics) quad-graph can be obtained 
by restriction to a subset of quadrilaterals in (cf. [7]), so this multidimensional solution 
yields a natural singular-boundary reduction of the type-Q ABS equations in that setting. Of 
particular relevance here is the further restriction to the odd (or even) sub-lattice which obtains 
the related Toda-type system on graphs; the singular boundaries here give the open boundary 
conditions for that system. 



13 Concluding remarks 

In this paper singular-boundary reductions have been defined for type-Q ABS equations and the 
most regular example, the singularity-bounded strip, has been studied in detail leading to an 
exact expression for its solution. The method exploits singularity structure in multidimensions 
to connect reduced systems with progressively fewer degrees of freedom. The basic idea is to 
reverse the chain, constructing the solution through an epitaxial sequence. Deciding whether or 
not this technique can be applied to other (perhaps less regular) singular-boundary problems 
comes down to a consideration of the geometric singularity configuration conditions (Theorem 
[T]) to see if such a chain can be constructed. 

The integration part of the problem has been solved by linearisation in multidimensions and is 
based on the procedure of obtaining soliton solutions developed in [33]. An important refinement 
needed here to incorporate the singularity structure is identification of the associated r function 
within that framework. The similar methods of construction allow to compare the solutions 
obtained here with the known soliton solutions. An important feature of the Q4 solitons is 
that their parameterisation involves elliptic and theta functions associated with a deformation 
of the elliptic curve of the equation. The solutions here are in this sense more intrinsic because 
they are in terms of elliptic and theta functions associated with the same underlying curve as 
the equation. On the other hand, these solutions are more complicated geometrically because, 
unlike the solitons, they are not defined isotropically throughout the lattice. In the expression 
for the solutions this manifests as a characteristic segmentation of the Hirota-type polynomial 
across the width of the strip. 
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